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Abstract 

We make a complete analysis of the Yukawa coupling unification in SU (5) extensions of the 
MSSM in the framework of the radiative symmetry breaking scenario. Both logarithmic 
and finite threshold corrections of sparticles have been included in the determination of the 
gauge and Yukawa couplings at Mz- The effect of the heavy masses of each model in the 
renormalization group equations is also included. We find that in the minimal SU(5) model 
b — r Yukawa unification can be achieved for too large a value of a s . On the other hand 
the Peccei-Quinn version of the Missing Doublet model, with the effect of the right handed 
neutrino also included, exhibits b — r unification in excellent agreement with all low energy 
experimental data. Unification of all Yukawa couplings is also discussed. 



March 1997 



E-mails : adedes@cc.uoi.gr, tamvakis@cc.uoi.gr 



I. INTRODUCTION 



Softly broken supersymmetry Jj]], possibly resulting from an underlying superstring 
framework, can lead to SU{2)i J x Z7(l)y g au ge symmetry breaking through radiative correc- 
tions for a certain range of values of the existing parameters || | |. A most appealing 
feature of the Minimal Supersymmetric Standard Model (MSSM) and its extensions is that 
it exhibits gauge coupling unification || in good agreement with the low energy values of 
the three gauge couplings, known to the present experimental accuracy, within the bounds 
set by the stability of the proton. The Higgs boson running mass-squared matrix, although 
positive definite at very high energies, when radiatively corrected develops at low energies 
a negative mass-squared eigenvalue that triggers electroweak symmetry breaking. Running 
parameters are studied in the framework of the Renormalization Group equations. 

In the present article we complete the analysis of electroweak symmetry breaking in 
the MSSM and various SU(5) extensions of it |J 0, focusing on fermion masses. Yukawa 
coupling relations, such as the b — r Yukawa coupling equality at high energies in relation to 
the experimental ^ value, give an extra constraint on models in addition to existing strong 
constraints like the experimental value of a s {Mz\ The input parameters of our analysis, 
apart from the standard low energy inputs (oem, Gf, Mz, quark and lepton masses) are 
the soft breaking parameters (M , Mi/ 2 , A Q ), the ratio of the two v.e.vs tan{3(Mz), the sign 
of the Higgs mixing parameter fx and the superheavy particle masses {Mh c , Ms,...) entering 
through their thresholds in the Renormalization Group equations. On the other hand, our 
output includes the strong coupling a s (Mz), the complete sparticle spectrum as well as other 
quantities, like the unification scale Mqut, etc. In this analysis, for obvious reasons, we shall 
concentrate on that part of parameter space that leads to a s compatible with experiment. 
The experimental values of third generation fermion masses will be introduced as input while 
the high energy values of Yukawa couplings will be treated as an output. 

We have included both logarithmic and finite one loop corrections (see Ref. |7j] for more de- 
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tails) in the calculation of sin 2 6{M z )\jy^ that determines the and ^Idr § au g e coupling 
values which are taken as boundary conditions at M z . Yukawa couplings are determined 
with the same next to leading order accuracy as we employed in our analysis for the gauge 



couplings. We extract the DR values of the Yukawa couplings at M z from the corresponding 
pole masses of quarks and leptons [§] as we describe below. 

The fermion masses, defined as the poles of the corresponding propagators, are related 



to the DR masses, rhf(Q), by the self energies, as 

m f le = rh f (Q) + $leX f (m p f de ) (1) 

where —iT,f(m P j° le ) = — + m^°' e S 7 ) is the one loop self energy on shell of quarks or 
leptons. Our results for self energies agree with those of references || JIIJ. In the case of 



the Peccei- Quinn version of the Missing Doublet Model (MDM+PQ) the presence of the right 
handed neutrino of mass Mr gives rise to a mass for the left handed neutrino m u = r " T ^ m/3 
via the see-saw mechanism [II]. The calculated one loop neutrino self energy correction is 



denoted by E^. From equation (1) we arrive at the values of Yukawa couplings at Mz 

Y * M *»™ ~ v(M z )\^cosf3(Mz) (2) 

(3) 
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where all hats mean running couplings or masses evaluated in the DR scheme. The running 
vacuum expectation value is given by the following formula, 

Mf + M| z (Mf ) = \(gl(M z ) + g\{Mz))v\Mz) (6) 



where M§ refers to the experimental Z-boson mass and H zz is the real transverse part of 
the Z-boson one loop corrections at Mz- 

For a given set of pole masses mt° le , ml° le , m p ° le we define the DR Yukawa couplings at 
Mz- Then we use the 2-loop Renormalization Group equations to run up to the scale Mgut, 
where g\ and g<i meet. We impose the unification condition 

9gut = gi(M GUT ) = g 2 (M GUT ) = g 3 (M GUT ) (7) 

and run down to M z . The thresholds of the heavy masses, both in gauge and Yukawa cou- 
plings, are treated with by using the well known method of the step function approximation 
[0] HI- We have assumed universal boundary conditions for the soft breaking parameters A Q , 
M Q , Mi/2- The whole procedure is iterated until convergence is reached, imposing the con- 
straints of radiative symmetry breaking, the proton decay bound, the experimental bounds 
on supersymmetric particles and the perturbativity of the Yukawa and gauge couplings up 
to M P . 

The value of the bottom quark pole mass can be determined indirectly through its influ- 
ence on hadron properties. In Table I we display the pole mass of the bottom quark as it 
is extracted from experiment in different processes following various techniques ||13|| . (For 
the determination of the pole mass of the bottom quark from QCD momentum sum rules 
for the T system see the recent article of Ref. |]14|.) 

TABLE I 

Description rrib(GeV) 
Lattice computation of T spectroscopy 5 ± 0.2 
e + e~ — > b hadrons cross sections 4.827 ± 0.007 
Heavy quark effective theory 4.61 ± 0.05 

T system 4.60 ± 0.02 

Table I: Values of the bottom quark pole mass 
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As a characteristic mean value of the bottom quark mass we shall choose below m& = 
4:.9GeV. The experimental mass of r-lepton is m T = 1.777GeV. In addition, we adopt 
for the mass of the top quark the average experimental value of .D0 and CDF experiment 
nit — 180GeU. We also take the r-neutrino to be in the cosmologically interesting domain 
m VT = 3 — lOeU. When we depart from the above values we indicate so. 

The unification of the Yukawa couplings in the MSSM or in some supersymmetric Grand 
Unified Models has been studied in || jn| (JTjJ. The presence of right-handed neutrinos as 



the only thresholds in the grand desert is studied in JT7J while some of their consequences 



in a supersymmetric 50(10) model in [IS]. In this article we examine, in a higher level 
accuracy framework, the ratios of the Yukawa couplings at the Unification scale Mqut m 
the minimal SU(5) and the Peccei-Quinn version of the Missing Doublet Model, in which 
the right-handed neutrino is included. 

II. A LITTLE BIT MORE ABOUT THE MINIMAL SU(5) MODEL 

In this section we apply the preceeding formalism on the minimal supersymmetric 577(5) 
model. We start with the superpotential |TP[ , 



W = -MxTriT?) + -^TriY?) + M 2 HH + X 2 HEH 
2 3 

+ V^Y^^H - \y^jH i,j = 1..3 (8) 

where the superfields E, H, H, 0, \1> belong to the representations 24, 5, 5, 5, and 10 
respectively. The SU(5) breaking occurs when the £ superfield develops a vacuum expec- 
tation value in the direction < E >= V Diag(2, 2, 2, — 3, — 3). In the resulting effective 
SU(3) x SU(2) x £7(1) theory we get the following heavy particles with quantum numbers: 

mS' 1 '"^ = Ml' 1 '^ = 5A 2 U (9) 

c -ti c 

Mq 8 ' 1,0 - 1 = M^ 1,3 ' 0) = M s = 5AiU , M^ lfi) = ^ (10) 



M v = 5g 5 V (11) 

In addition to the RGE's for the gauge couplings 0, we display here the RGE's for the 
Yukawa couplings, treating the superheavy thresholds in the step function approximation 

16vr 2 ^ = Y b {6Y 2 + Y 2 + Y T 2 + 3Y 2 Q(H C ) + 2Y 2 Q{H C ) + ^A 2 2 9(u,) 

+ ^XlQ(S) + 3X 2 2 Q(V) - (Lfi + 3g 2 2 + fg 2 3 ) - 8g 2 5 Q(V)} (12) 



= Y T {4Y 2 + 3Y 2 + 3 -X 2 2 Q(u) + ^X 2 2 Q(S) + 3X 2 2 Q(V) 

+ 3Y 2 Q(H C ) + 3Y 2 Q{H C ) - ( 9 -g 2 + 3g 2 2 ) - 12g 2 5 e(V)} (13) 



AY Q 

16tt 2 ^ = Y t {QY 2 + Y 2 + 3Y 2 Q(H C ) + 3Y 2 Q(H C ) + -X 2 2 &(S) 

+ lx 2 2 Q(uj) + 3X 2 2 &(V) - (^g 2 + 3g\ + fg 2 3 ) - 10g 2 Q(V)} (14) 



16vr 2 ^ = 3X 1 Q(u;){^X 2 i e(S) + 3A 2 0(V) + A 2 - Ag 2 2 - QglQ{V)} (15) 



16tt 2 ^ = X 2 Q(uj){3Y 2 + Y 2 + 3K t 2 + 3A 2 0(cu) + ^X 2 2 0(S) + 6A 2 0(V) 

+ ^A 2 6(S) + 3X\Q{V) + A 2 - (^ 2 + 7g 2 2 ) - 12<7 5 2 e(V)} (16) 

where t = Inj^. We have denoted Q(x) = 0(Q 2 - M 2 ). Note that Mqjjt = 
max{My, Mh c , Af s } which we take equal to My as it turns out to be in all relevant cases 



In Figure 1 we have plotted the ratio of bottom and tau Yukawa couplings at Mqut as a 
function oitanfi and the input high energy mass Ms. We select as input soft breaking masses 
at MquT) so as to give consistent low energy results within the experimentally acceptable 



region. We obtain that the b — t unification in the minimal SU(5) is in agreement with the 
pole mass of the bottom quark extracted from the Lattice computation of T spectroscopy as 



it is shown in Table I. Nevertheless, even in this case the output MS value of a s (Mz) is quite 
large (a s (Mz) > 0.127) compared with its average experimental value 0.118 ± 0.003 |L3| . 
Both the proton decay bound and the theoretical perturbativity bound on the top Yukawa 
coupling (Y t < 1.5 at Mqut) strongly constrain the Yukawa unification. We have also seen 
that Yukawa unification is sensitive through a s to the logarithmic and finite self-energy 
corrections of Yukawa couplings as well as to the analogous corrections to the Weinberg 
angle at Mz- As it is obvious from Figure 1, decreasing the value of Ms destroys (for the 
lower values of m&) or restores (for higher values of mj) the picture of Yukawa unification in 
contrast to the MSSM case, in which all the high energy particles have been decoupled from 
the renormalization group equations. In Ref. the hope was expressed that lowering the 
value of Ms might result in the unification of gauge couplings at a scale close to the so-called 
string scale M str i ng = 5 x 10 17 GeV. This can be the case only if mj > 5GeV, as shown in 
Figure 1, but again with too large a value for a s 0. A proposed mechanism to reconcile 
the value of a s with the experimental data is the one which goes through the high energy 
threshold of colour triplet mass Mh c - Lowering this value down to \Q lb GeV we obtain values 
of a s (Mz) which are in good agreement with the experimental ones. This cannot be the 
case however in minimal SU(5) where we encounter the bound M# c > 1.8 x 10 16 GeV due 
to the proton decay constraint. This problem is solved in the Peccei-Quinn version of the 
Missing Doublet Model where the proton decay rate is suppressed due to the presence of the 
intermediate scale 10 12 GeV . Putting aside the problem of large a s in the minimal SU(5) 
model and varying M# c in the region (1.8 — 3) x 10 16 GeV we obtain values for the b — t 
Yukawa ratio, inside the shaded region of Figure 1. 
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III. THE PECCEI-QUINN VERSION OF THE MISSING DOUBLET MODEL 



In order to avoid the numerical fine-tuning required in the minimal 577(5) model we can 
replace the adjoint Higgs representation by the 75 representation which couples the Higgs 
pentaplets to an extra pair of Higgses G and Q in the 50+ 50 representation respectively 
that contains no isodoublets. This extension of the minimal SU (5) is known as the Missing 
Doublet model |2l| . The superpotential of this model is 

W = MxTr(S 2 ) + iA!Tr(S 3 ) + A 2 /7£0 + A 2 i7£0 + M 2 00 

o 

- lY&V&iH + ^VffirfF U = 1-3 (17) 

As it has been shown in Ref. J7j this model built in such a way, although it predicts small 
values of a s (Mz), leaves a narrow window in the parametric space (especially in the tanfl) to 
work with, because of the bound which comes from the proton decay stability. This problem 
which is also present in the minimal S77(5), albeit with large values of a s (Mz), provided 
strong motivation to construct versions of S77(5) with a Peccei-Quinn symmetry [EI3] that 
naturally suppresses D = 5 operators by a factor proportional to the ratio of the Peccei- 
Quinn breaking scale over the GUT scale [^J. If we promote the symmetry by an extra 
global Peccei-Quinn U(l) factor the terms responsible for Higgs masses must be replaced by 

X 2 Hm + A 2 77£0 + X' 2 H'm' + x' 2 H'm' 

+ M 2 ee' + m 2 q'q + x 3 ph'h' 

We have introduced a new set of chiral supernelds,if'(50), if'(50), 6'(50), (50) and P is a 
gauge singlet superfield. The charges under the Peccei-Quinn symmetry are 0(— 1/2), 

H(-2), #(-1/2), 0(2), 0(1/2), 0'(-2), ©'(-1/2), T?{l/2), H'(2), P(-5). The Peccei- 
Quinn symmetry is broken Q when the gauge singlet P develops a vacuum expectation value 



"This is achieved through a suitable sector involving singlets [24]. 



7 



at an intermediate energy < P >= ~ 10 10 — 10 l2 GeV. The model is naturally extended 
by introducing the three right-handed neutrino superfields N? adding to the superpotential 
the following Yukawa interactions: 

Wiv = Y u N c (j)H + ^\ 4 N C N C P (18) 

These terms induce Majorana masses for the right-handed neutrino multiplets iV? and these 
masses are proportional to the intermediate breaking scale since Mr = ^Mh> ~ O(10 n — 
10 12 )GeV. Finally, the terms incorporated in Wat induce through the see-saw mechanism 
[|Tl"|1 very small masses for the neutrinos. The spectrum obtained after the breaking of the 
SU(5) gauge group can be read from Table II where we have used the following reasonable 
assumptions for the parameters in the superpotential of this model: 

A 2 ~ A 2 ~ A' 2 ~ A 2 M 2 ~ M' 2 ~ 10 18 Ge1/ (19) 



TABLE II 


75 


50 


Colour Triplets 


M 58 (8,3,0) = M E 


M(6,l,|) =4M 2 


M Hc ~ Mflj ~ 32Aig 


M 5 (1,1,0) = 0.4M S 


M(8,2,|) =8M 2 


M 6c ~ M e , c ~ 8M 2 


M 56 (6,2,|) = iV% 6 (6,2,-|) = 0.4M s 


M(6,3,-|) =4M 2 


Doublets 


M o (8,l,0) = 0.2M S 


M(3,2,-|) = 8M 2 


M Hf =0 


M 5a (3, 1, |) = A% 3 (3, 1, -|) = 0.8M S 


M(l,l,-2) = 8M 2 


M H > = A 3 < P >~ 10 n GeV 


M E = 20Mi = -f Ail/ 




My = V6g 5 V 



Table II: The spectrum of the MDM+PQ model 



The running of the Yukawa couplings at the one loop level is given by the following 
renormalization group equations, 
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16vr 2 ^ = Y b {6Y b 2 + Y? + if + \y 2 Q{H c ) + \y?Q{H$ + Y?Q{H e ) 

7 1 fi 

+ Y 2 Q(K) + Y U 2 Q(H C )Q(R) - {-g\ + 3gl + -fgl) - 8g 2 5 Q(V)} (20) 



16vr 2 ^ = Y T {4Y T 2 + 3lf + ^ 2 6(tf c ) + \yMK) + \^{H C ) 

+ \yMK) + YMR) ~ (~<7i 2 + ^l) - 12^ 2 6(V)} (21) 



+ \y?Q{H> c ) + Y 2 Q{R) ~ (^i 2 + H + fgl) - m?MV)} (22) 



= Y u e(R){Y T 2 + 31? + 2Y 2 + 3 ~Y*e(H c ) + \y 2 Q{H' c ) 

+ 2i; 2 e( J R) + 3i; 2 e(F c ) - (^ 2 + 3^ 2 )} (23) 



16vr 2 ^ = 3A 1 0(O){^A? + ^A 2 6(S) + 144A 2 6(S 6 ) + ^A 2 6(S 3 ) 

+ ^A 2 e(s 8 ) + ^A 2 e(v) - 6^ 2 - ios 5 2 e(\/)} (24) 

Note that O(R) stands for Q(Q 2 - M|) and that again we consider cases with My = Mgut- 

We have taken for all figures arranged below the input soft breaking masses at Mqut to 
be: A Q = 400GeV, M = 300GeV and M 1/2 = 300GeV\ Also the superheavy mass M E is 
taken to be M E = 10 16 GeV\ These values are indicative values corresponding to acceptable 
a s 0. 

In Figure 2 we have plotted the ratio of the b — r Yukawa couplings at Mgut as a function 
of tanf3(Mz) when we vary the pole mass of the bottom quark in the range allowed by Table I. 
It is clear in the case of the MDM+PQ model we get b — r Yukawa unification both for \x > 
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or /i < 0[| together with the excellent agreement of the value of a s [Mz) as this is compared 
with the experimental average value 0.118 ±0.003. In the case of fi > for large tan/5 there 
are significant finite threshold corrections in the bottom and tau Yukawa couplings which 
arise from the chargino/squark or gluino/squark loops and they are proportional to [xtanf}. 



This fact [[L5| supports b — r Yukawa unification in the large value region of tan(3 for a > 0. 
Note also that the upper bound in the input parameter tanft(Mz) < 36 does not come from 
the bound of the proton decay rate as in the case of the minimal SU(5) model, but from 
the requirement of the radiative symmetry breaking. This is a special feature of the Peccei- 
Quinn version of the Missing Doublet Model and is valid for all figures arranged here. 

In Figure 3 we display the ratio of both b — r and t — b Yukawa couplings as the top pole 
mass is varied inside the region mS° t : = 182 ± 8GeV in the MDM+PQ model. The £ ratio 

it 

then takes values in the regions 0.91 — 1.09 for a < and 1.00 — 1.11 for a > respectively. 
Note also that as can be seen in Fig. 3c and Fig.3d the t — b ratio never reaches unity due to 
the constraint of the radiative symmetry breaking. 

Proceeding next to Figure 4 we vary the mass of the colour triplet M# c (Fig4a,b) and 
the mass of the intermediate doublet Mn' f (Fig4c,d) for a < and > keeping fixed 
the mf ie = ISOGeV and m p b ole = 4.9GeV. Acceptable values for a s are obtainable if we 
stay in the vicinity of M# c = lQ 15 GeV and Mr* = 10 ll GeV 0. For larger values of the 
colour triplet or smaller values of the intermediate doublet we get larger values for the strong 
coupling a s completely ruled out by experiment. As can be seen in Figure 4 6 — r unification 
seems to prefer the case where /i < and low values of tan/3 ~ 5. If we lower the ml° le down 
to 4.8GeV we also get for /i > (this is shown in Fig. 2) b — r unification in the low and 
large tanf3 regime. 



' We follow the conventions of Ref. 
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In Figure 5 we consider the behavior of all ratios of Yukawa couplings. The Renormal- 
ization group equations for the three ratios Rt- V = y-, Rb-r = y 1 and R t -b = ^ are 
dRt- v 



dt 



= R t -„{3(Y t 2 - YMR)) + W " *?) + 3(^ 2 - Y U 2 )Q(H C ) 
A 16 

;^ + y^ 2 )-lO^ 2 0(^)} (25) 



dR b - 



T = J R fe -.{3(r b 2 - r r 2 ) + (if - yMR)) + (*? - ^ 2 )0TO 

:y£ 3 2 -^ 2 )+% 2 e(V0} (26) 



- RtMHY? - n 2 ) - tt? - n 2 ©(^)) + (^ 2 - y^m 

- 2 -gf - 2g!&(V)} (27) 

It can be easily observed that the 5'O(10)-type condition Y t = Yb = Y T = Y v in the limit 
of high energies where all thresholds theta functions are equal to 1 corresponds a zero of 
the beta functions for the ratios provided we neglect the contribution of gauge couplings. In 
contrast to that the SU (5) condition Y b = Y T leads to a zero of the /3r 6 _ t = dR ^ T function 
in the presence of gauge couplings. Note that neither of the above occurs in MSSM where 
the contributions of the extra fields imposed by the extended GUT symmetry are absent. 

Figure 5 shows that the S'O(10)-like Yukawa unification requires as expected large values 
of tanf3 which are not allowed by the requirement of radiative symmetry breaking. Note 
however that this is not necessarily true in a realistic 5*0(10) model. Although complete 
Yukawa unification is excluded YJ, = Y T and Y t = Y u can be realized independently for a wide 
range of tanf3 values. It should be noted that the neutrino pole- mass that corresponds to t — v 
Yukawa coupling unification is approximately m v °^ e = (lOOGeV) 2 / Mr. This, very roughly, 
corresponds to a low energy value of the top Yukawa coupling approximately twice that of 
the neutrino Yukawa coupling. For a cosmologically interesting neutrino of m v °^ e = lOeV 
this fixes the lepton number-PQ-breaking scale to be Mr = 10 12 GeV. 
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dt 



IV. BRIEF CONCLUSIONS 



In this article we have studied in a 2-loop accuracy framework the unification of the 
Yukawa couplings in extensions of the MSSM which are based on the group SU(5). Both 
logarithmic and finite threshold contributions of SUSY particles have been included in the 
calculation of the gauge and Yukawa couplings at M z . We have taken into account the effect 
of all heavy masses in the renormalization group for the gauge and Yukawa couplings. The 
models that we have examined are the minimal SU(5) model and the Peccei-Quinn version 
of the Missing Doublet model. In the minimal SU (5) model b — r unification is possible 
for m b values larger than 5GeV in the low tanf3 regime but in all cases the encountered 
a s (Mz) values are too large to reconcile with experiment. Due to the bound of the proton 
decay stability we cannot go to the large tanf3 regime in this model. On the other hand, 
in the Peccei-Quinn version of the Missing Doublet model, where the effect of the right 
handed neutrino is naturally included, we have found (Fig.2,3,4) that b — r unification at 
Mqut can be achieved both in the low (/i > and \x < 0) and the large (// > 0) tan(3 
regime for a s values compatible with the low energy experimental data. It is interesting to 
note that this is achieved for value of the intermediate Peccei-Quinn b — r lepton number 
breaking scale in the cosmologically interesting neighborhood of 10 12 GeV. We have also 
considered whether all the Yukawa couplings are unified at the scale Mqut (Fig.5). Only 
b — t Yukawa coupling unification is possible within this model since the constraint of the 
radiative symmetry breaking forbids an SO (10) like condition Y t = Y b = Y T = Y u at M GUT . 
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Minimal SU(5) 




tan|3 



Figure 1: The allowed region for the ratio of bottom and tau Yukawa couplings at Mqut- The input 
values of the tau lepton and top quark are m T = X.lllGeV and m t — 180GeV respectively. The extracted 
value of a s (Mz) is greater than .127 . If we still lower the My, value down to 10 14 GeV, we observe quite large 
values of a s . In the upper curve of the shaded regions the heavy particles Ms, Mh c have been decoupled 
(MSSM curve). 
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MDM+PQ Model 

1.30 i — i — i — i — i — i — i — i — i — i — i — i — 




Q YQ I I I I I I I I I I I I I I I I I I I I I 

4 8 12 16 20 24 28 32 36 40 

tan (3 

Figure 2: The ratio of the b — r Yukawa couplings as a function of tan j3 in the MDM+PQ model when 
the m^ ole is varied for fi < and /i > 0. The fixed input parameters are mP t ole — l80GeV, M H ' f — lO n GeV, 
M Hc = 10 15 GeV, M R , = W n GeV, m? v ° le = 10eV. 
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Figure 3: The ratio of the b — r and t — b Yukawa couplings at Mqut as the top pole mass is varied. We 
keep the bottom pole mass, mP^ le = 4.9GeV fixed. The horizontal line shows the b — r or t — b unification. 
The other inputs are similar to those displayed in the Figure 2. 
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Figure 4: The effect of the variation of the Mh c and Mh' , within the range allowed by proton stability 
in the b — r unification and also in the extracted value of a s . We take as inputs: m£°' e = 4.9GeV and 
m p t ole = 180GeV. 
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Figure 5: Unification of all Yukawa couplings in the MDM+PQ model with m p b ° le = A.9GeV, m\ ole = 
180GeV and the other input parameters displayed in Fig. 2 
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